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In recent years quantum Monte Carlo methods have become powerful numerical tech- 
niques for attacking the many-body problem in condensed matter, nuclear physics, and 
other fields. Unfortunately, these methods are also suceptible to what is known as the sign 
problem. Generically speaking, the sign problem arises when an expectation value to be com- 
puted is the average of a wildly fluctuating quantity, so that numerical noise overwhelms 
the signal. 

This paper focuses on the auxiliary-field path integral, which is evaluated via Monte Carlo 
sampling, and further concentrates on the application to the nuclear shell model ||l]-§. The 
nuclear shell model is an excellent venue for auxiliary-field Monte Carlo (AFMC) because 
one can prove that a large class of systems are free from the sign problem. The 'goodness' 
of an interaction with respect to the sign problem is related to the time-reversal properties 
of the effective one-body Hamiltonian |l[]. 

This paper derives two sum rules for the time-reversal properties, and thus the sign- 
problem properties, of shell-model interactions. The sum rules are ultimately expressed in 
terms of the pairing matrix elements and prove two points of conventional, but previously 
unproved, wisdom in AFMC computations. 

I. AUXILIARY-FIELD PATH INTEGRALS 

Traditional shell-model calculations diagonalize a Hamiltonian if in a many-body basis; 
the current best limit on model-space dimension is about 5 x 10^ [Q although more typical 
spaces have basis dimensions of 10^-10^. Rather than diagonalization, AFMC shell model 
calculations employ the imaginary-time evolution operator exp{—f3H) to project out the 
ground state and evaluate thermally weighted expectation values. For the purposes of this 
paper is it sufficient to sketch the broad outlines of the method; for a detailed description 
the interested reader is referred to |[1|,0]. 

The nuclear Hamiltonian H is usually two-body and can be written as quadratic in 
one-body operators 



H = Y.,^p^ + Wyo.pl. (1) 

This is done in detail below. The quadratic term is source of all troubles; in order to facilitate 
numerical computation it is linearized, via the Hubbard-Stratonovich transformation, to 
form an effective one-body Hamiltonian that depends on the imaginary time r: 

^(^) = Xl(^a + SaVa(Ja(r))pa, (2) 

a. 

where Sq = ±1 if V^ < and = ±i if V^ > 0. Then the imaginary-time evolution operator 
exp(— /3i^) becomes a path-integral: 



exp(-/5^)=y'2)[a]expU^^''dr^|K.|V„(r)2j T exp (- J^^ drhir 



(3) 



where T denotes time-ordering. This path integral can be evaluated by Monte Carlo sam- 
pling if a sufficiently well-behaved weight function W{a) can be found. An ideal weight 
function could be, for example, the trace of the integrand of Eqn. (|^) itself. Unfortunately 
this choice of W{a) often takes on both positive and negative values and one must use \W\ 
as a weight function. In those cases typically (W^/|W^|) <^ 1, and the statistical errors for 
the Monte Carlo value of the observable are overwhelmingly large. For AFMC methods this 
is the manifestation of the sign problem. 

II. APPLICATION TO THE SPHERICAL SHELL MODEL 

The shell model is built upon single-fermion states, using fermion creation and destruc- 
tion operators a^, a. For a spherical basis these fermion states have good angular momentum 
j and m, as well as orbital angular momentum i which dictates the parity (—1)^. For each 
state we associate a "time-reversed" partner, j and — m, and define a time-reversal opera- 
tion by a bar: dj^^ = (— l)"'^'"^^ctj,-m and aj ,^ = (— l)-'+™+^a] _^. Because of the half-spin 
statistics, dj^n = —a-j,m- As hinted above, the time-reversal properties play a central role. 

The two-body Hamiltonian is usually written as 



^ abed JT MT, 



(4) 



where (at = Vl + ^ab, the Vjt are matrix elements of the interaction, and the two-body 



creation operator is Aj 



JT,MT^ 



(ab) 



Ja Jb 



. reforming a Pandya transformation 



JT,MT^ 

from the particle-particle to the particle-hole representation, the Hamiltonian becomes 

-f^2 = - XUZ ^Kiiab, cd) Y^{-1)^ PKM,Mc)pK-M,i{hd), (5) 

^ abed KI M 

plus a residual one-body term that does not concern us here. The one-body or density 



operator is Pkm{0'C) 



«l ® «> 



KM 



. and proton and neutron densities are combined: 



PKM,I = PKM,p + ( — 1) PKM,n- 



(6) 



Note: We will refer to these densities as 'isoscalar' (/ = 0) and 'isovector' (/ = 1) even 
though the latter do not truly have good isospin. 

The matrix elements in the particle-hole representation are 

ja jb J 



iy^+HabCcdT.i-^y i'^J + ^) 

J I jd jc K 



EK,i=o{ac, hd) = 

X 2 {KT=ii(^b, cd) + - [Vj%^o{ab, cd) - Vj%^^{ab, cd) 

ja jb J 



(7) 



Ek,i=i{0'C, hd) 



J lid jc K 

X - \y^T=o{ah, cd) - VfT=iiab, cd) 



(8) 



V is a general two-body matrix element which can be separated into symmetric and an- 
tisymmetric parts, V'^ = V^ + V"^; the symmetry properties are given by Vjj- (ab, cd) = 
^[Vj^{ab,cd) ± (— l)"^+-5'"+-'''+-^~^Vj^(6a, C(i)]. Of course, only the antisymmetric matrix el- 
ements V^ are physical and are used in Eqn. (^; the symmetric matrix elements V^, 
although unphysical, are a degree of freedom in the particle-hole representation and can be 
used to encourage convergence. 

The next step is to diagonalize the Eki, which has eigenvalues Xk,i{c() and eigenvectors 
VKi,a{o,c) ■ Define Pki{c() = Y.ac'^Ki,a{cic)pKi{(ic)] then finally 



■^ KI a M 

Now the Hamiltonian is in the form of a sum of manifestly quadratic operators and the 
Hubbard-Stratonovich transformation can be apphed. 



III. LANG'S THEOREM AND SIGN-PROBLEM-FREE INTERACTIONS 

In some systems W is positive-definite and there is no sign problem. Lang's theorem |lj] 
shows that for nuclear physics a wide class of systems are free from the sign problem. 

Any one-body operator p can be separated into time-even and time-odd parts: p = 
^{p+p) + ^{p—p). The basic statement of Lang's theorem is that, for even A^, even Z systems 
(or N = Z odd), if the effective one-body Hamiltonian (^ is strictly time-even, then there 
is no sign problem; specifically, that (\E'r| Texp (— j/j dTh{a)) |\E'j') > 0. Since i = y/—T is 
also time-odd, this means that time-even operators pa must have a real coefficient and time- 
odd operators an imaginary coefficient. Because imaginary coefficients only enter ii Va > 0, 
this means that time-even terms arise from attractive interactions and time-odd terms from 
repulsive interactions. One can easily show that Pkm{(ic) = (— l)'^'^*'^7r(ac)px,-M(a.c), where 
7r(ac) = parity of density operator = (—1)^"+^'=. (The parity can only be negative when 
interactions cross major harmonic oscillator shells f^. For a Ohu model space, such as the 
sd or pf shells, the parity can be ignored.) Hence the criterion for a 'good' interaction is 

{-l)^n{a)XKiia) < 0. (10) 

In order to deal with realistic interactions 0,^ which modestly violate Lang's rule, Al- 
hassid et al. introduced a method whereby one deforms the realistic interaction to one 
that conforms to Lang's rule and then extrapolates from the "good" regime to the "bad" 
regime. That is, let 

H = i/good + -f^bad (11) 



where Hgood satisfies Lang's rule and -ffbad uniformly violates it. This is done by decomposing 
the interaction into the form (^), and segregrating 'good' from 'bad' interactions via (|10|). 
Then one introduces 

H{g) = f{g)H^ood + gH^^d, (12) 

where one usually takes f{g) = [1 — (1 — 5')/4]. For g = 1 we have the original Hamiltonian, 
but for g < 0, H{g) satisfies Lang's rule and is amenable to the Monte Carlo shell model. 
The typical procedure is to compute for several values oi g < and then extrapolate to 
g = 1. The results agree for model spaces where exact diagonalizations can be performed 
for comparison p|,||. 

There are two points of conventional wisdom regarding Monte Carlo Shell Model calcu- 
lations, both of which are noted in Lang et al 0] and later papers |@J^: 

Conventional wisdom ^ 1: It is customary to set the unphysical Vfrp^^ to Vfj^^Q so 
that Ek,i=i = 0. This reduces the number of auxiliary fields to be integrated over. On the 
other hand, it is conceivable that for some interactions one could exploit Ek,i=i 7^ to avoid 
the sign problem at the modest cost of doubling the number of auxiliary fields. 

This paper proves explicitly that the latter statement is impossible and that the con- 
ventional choice is indeed the best, not only with regards to reducing the dimensions of the 
integral but with regard to alleviating the sign problem. 

Conventional wisdom ^ 2. The extrapolation procedure amounts to increasing the 
strength of the pairing interaction. Ref. ^ notes that an attractive pairing interaction, 



with V}io 5.^1(00, cc) = -G^{2ja + l)i2jc + l) yields Xk,i=o = -{-l)^G/2 for all K; if 
added crudely in sufficient strength to any interaction it will make it 'good.' Empirically 
one finds that the more general extrapolation procedure described above seems to increase 
the pairing strength @J^. The sum rules given below show explicitly a relation between the 
sign problem and the pairing interaction. 



IV. TWO SUM RULES AND THEIR IMPLICATIONS 

This is the central section of the entire paper. The sum rule is, for isoscalar densities 
(/ = 0) is 

Y,i2K+l)J2{-lfnia)XK,i=oia) = W ^(«c) V(2j. + l)(2j, + l)Vjto,T=iiaa, cc), (13) 

K ce ■^ ac 

and for isovector densities(/ = 1) 

J2{2K + 1) J2{-lfn{a)\K,i=i{a) = 0. (14) 

K a 

The proof of these sum rules is straightforward. For a given K, /, Y.a ^ki{.o) = 
J2acEKi{ac,ac), since the sum of the eigenvalues is just the trace of the matrix. Then, 
using 



ja jc K I (^_l-^K+ja+jc 



3c Ja J ^(2j,+ l)(2j, + l) 
and the orthogonality of six-j symbols [0 , one finds 



j:i2K+l)i~lfY.>^K,i=o 



a) 



K 



2 Y. (K7lo,r=i(aa, cc) + 2 \^Lo,t=q{<^<^. cc) - 1/jlo,r=i(aa, cc)] j , (15) 

^(2ir + 1)(-1)^ Y. >^K,i=M) = -\Y: [K7to,T=o(««, cc) - Vf=o,^=i(aa, cc)] . (16) 

K ct ^ ac 

(For each of these steps parity must be conserved so that implicitly 7r(a) = 7r(ac) in equations 
(p!5|), ([T6|). ) The final step is to note that, due to their symmetry properties, Vj=o,T(oa, cc) 
must be symmetric for T = and antisymmetric for T = 1. Hence the terms VjIq rp^^i^aa, cc), 
Vj^Q rp^i{aa, cc) vanish and VJLq rp^i{aa, cc) reduces to Vj^Qrp^i{aa,cc). Thus the final re- 
sults follow easily. 

Now to interpret these sum rules. Recall that Lang's rule dictates that a 'good' interac- 
tion will have (—1)^ n{a)\Ki{a) < 0. By inspecting the 'isovector' sum rule, if any Xk,i=i 
is 'good' (and nonzero) there must be another, nonzero 'bad' A to cancel its contribution. 
The only option to avoid the sign problem in the isovector channel is to set all Xk,i=i = 0, 
thus confirming conventional wisdom ^ 1. 



The isoscalar sum rule is expressed in terms of the J = pairing matrix elements. First 
consider 'normal parity' pairing matrix elements, that is, within a major oscillator shell. 
If Hbad —>■ —Hbad, then the contribution of the 'bad' A's will switch signs and increase the 
magnitude of the isoscalar sum rule. This must then manifest as stronger, more attractive 
pairing matrix elements, as expounded in conventional wisdom #2. The 'abnormal parity' 
pairing matrix elements, which takes a J = pair of nucleons from one oscillator shell to 
one of opposite parity, must however become more repulsive. 

Such behavior is illustrated in Figure 1 for a realistic cross-shell interaction |^. There 
we plot the J = matrix elements Vj=oT=i{ab, cd), grouped by shells, denoting the srf-shell 
by + and j9/-shell by — . The case g = 1 (filled circles) correspond to matrix elements from 
the original, realistic interaction that has a 'bad' sign, while g = (open circles) correspond 
to an interaction minimally deformed to have a 'good' sign. As described by the isoscalar 

sum rule, the intershell (V(++, ++) and V{ , ) ) pairing becomes more attractive as 

g = 1 —^ 0, while the intrashell (F(++, )) becomes more repulsive. Matrix elements 

of the type V{-\ — , H — ) do not exhibit any significant trend, which is accord with the sum 
rules: such matrix elements do not contribute to the trace over the Ex{ac, ac). 

In summary, the two sum rules derived herein allow one to confirm the conventional 
wisdom regarding AFMC computations. They also potentially yield a powerful tool to 
analyze the extrapolation and to look for alternatives. This is left to future work. 

This work was done with the support of the Department of Energy under grants DE- 
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by Lockheed Martin Energy Research Corp. for the U.S. Department of Energy under 
contract number DE-AC05-96OR22464. 
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FIGURES 
FIG. 1. J = pairing matrix elements for realistic interaction in sd-pf shell model calcu- 
lations, g = 1 corresponds to original, realistic Hamiltonian that has a 'bad' sign; (7 = is an 
interaction minimally deformed to have a 'good' sign, x-axis is arbitrary ordering of J = matrix 
elements. '+' denotes the sd shell while '-' denotes the pf shell. 
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